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Abstract 

The real part of the optical potential for the nucleon-nucleus scattering 
at lower energies (E{ <100MeV) has been calculated including nucleonic 
and mesonic form factors by a double folding approach. Realistic density- 
and energy-dependent effective iViV-interactions DDM3Y, BDM3Y and 
HLM3Y based on the Reid and Paris potentials are used in this respect. The 
effects of the nucleon density distribution and the average relative momentum 
on the folded potential have been analysed. A good agreement with the phe- 
nomenological potential of Lagrange-Lejeune, as well as with the parametriza- 
tion of Jeukenne-Lejeune-Mahaux for both neutron and proton double- folded 
potentials is obtained. The results indicate that the strongly simplified model 
interactions used in preequilibrium reaction theory neglect important dynam- 
ical details of such processes. 
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The nucleon-nucleus optical model potential (OMP) is an important tool for under- 
standing and interpreting the mechanism of the interaction between a nucleon and a nuclear 
system. The explicit forms of the t-matrices and nuclear matter g-matrices have been used 
through the folding procedures to obtain more realistic OMP for practical use in the nuclear 
reaction calculations |1|]. Brieva and Rook || used the local density approximation (LDA) 
in order to apply to finite nuclei the nucleon-nucleon (NN) g-matrix effective interaction 
appropriate for a pair of nucleons interacting in infinite nuclear matter. Jeukenne, Lejeune 
and Mahaux ( JLM ) || recommended a simpler parametrized form of the isoscalar real com- 
ponent obtained from Brueckner-Hartree-Fock nuclear matter calculations, to be used for 
practical nuclear reaction calculations. Actually, the use of the free nucleon-nucleon t-matrix 
as a complex effective interaction has already been applied successfully to nucleon scattering 
at higher energies (Ei >100MeV), e.g. ||] and references therein. 

Microscopic descriptions have become a standard approach for the description of re- 
actions leading to discrete final states and giant resonance, respectively [Q. However, a 
different situation is encountered for more complex nuclear reactions involving multistep 
processes. In this paper, we consider the inclusion of a microscopic OMP into the descrip- 
tion of preequilibrium emission processes. Hitherto, the complexity of such reactions has 
been considered to inhibit the use of microscopic methods. It is common practice to use 
strongly simplified models, for example the simplest 1 fm range Yukawa interaction has 
been chosen to describe the effective A^AMnteraction leading to particle-hole excitations 
in the frame of the Feshbach-Kerman-Koonin (FKK) model ||, its strength Vq being the 
only free parameter adjusted to reproduce the experimental data. The large discrepancies 
still present at low incident energies (Ei <50 MeV) in the systematics of the phenomeno- 
logical Vo values prove the necessity to consider a more realistic ATAMnteraction which 
should be chosen consistently with the corresponding real part of the optical model poten- 
tial. One may consider in this respect the density- and energy-dependent effective DDM3Y, 
BDM3Y and HLM3Y interactions based on the g-matrix elements of the Reid and Paris 
AW-potentials [|3|-|T0|j, which seem to be most successful for the nucleon-nucleus scattering 
within the low-density nuclear matter ||. 

The present work concentrates on a comparison of folding OMPs from different inter- 
actions and density models. The intention is to derive a global model which describes 
reliably elastic scattering over a large range of incident energies and target masses. Once 
such a model is available, it can be applied on safer theoretical grounds to the various re- 
action channels of preequilibrium reactions. Obviously, the OMP in these channels cannot 
be determined directly by experimental means. A microscopic model will reduce the still 
persisting uncertainties. Moreover, the consistency between the interaction used in elastic 
scattering and the non-elastic processes is guaranteed. 

Usually, effective NN interactions are parametrized in terms of meson exchange-type 
propagators, e.g. of Yukawa shape. Folding calculations for OMPs and transition potentials 
involve an averaging over the intrinsic momentum distribution of the target nucleons. Hence, 
a region of off-shell momenta is covered. Under such conditions the use of "on-shell" meson 
exchange propagators might lead to spurious effects because contributions from large off- 
shell momenta are overestimated. These effects are well-known from calculations of the free 



NN t matrix [13] and in Brueckner theory, e.g. 14,15]. In order to avoid these uncertainties 



one has to introduce vertex and nucleon form factors which regularize the large momentum 
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region. Their effect is to exclude contributions from small distances where the composite 
structure of mesons and nucleons would become visible. 

For a OBE-type interaction the momentum space vertex cut-offs correspond to remove 
in coordinate space the singularities of Yukawa functions at the origin. The same effect is 
obtained by folding an r-space form factor to a Yukawa function. In our calculations we 
follow this prescription and use a double folding procedure with effective NN interactions 
and a form factor describing the intrinsic structure of the nucleon. We expect to take into 



account typical features of realistic NN potentials (e.g. [pT|-[T3]|) of which the regularized 



Reid93 potential [I2| includes explicitly a dipole form factor. 

The main aim of this work is to study the effects of both nuclear and nucleon density 
distributions. Equally important is a reliable treatment of the non-local exchange contribu- 
tions from anti-symmetrization. In order to avoid the complications of a non-local optical 
potential we use the local momentum approximation. Different prescriptions of the related 
quantities like the kinetic-energy density r(r) and various approximations of the average 
relative momentum k av (r) [^6| , |17| , on the corresponding double- folding real OMP at lower 
incident energies (Ei <100 MeV) are investigated in detail. 

The general expressions for the direct and exchange double-folding (DF) real parts of 
the optical potential in terms of the nuclear densities pi(ri) and p 2 (r 2 ) of the projectile and 
the target nucleus, respectively, and the effective iViV-interaction v (s=R + r 2 - ri) are by 
assuming proportional proton and neutron densities |],|8| , 



^O(l) 



(E,R) = 5 (i)<5 (2) / dri J dr 2 pi(ri) p 2 (r 2 ) v^ 01) (p,E,s), (la) 



Uq*)(E,R) =<y (i)£o(2) / dr i J dr 2 pi(r 1 ,r 1 + s)p 2 (r 2 ,r 2 - s) VQ f 01) (p, E,s) exp 



ik(R)s' 

M 

(lb) 



where 5 =1, 5i=(Nx - Z 1 )/A 1 , 5 2 =(N 2 - Z 2 )/A 2 , M=A 1 A 2 /(A 1 + A 2 ), k(R) is the incident 
relative momentum, and pi and p 2 in Eq. ( flb"|) are the density matrices for the projectile and 
the target nucleus, respectively. If the projectile is a nucleon, in the point-like approximation, 
i.e. for p(r)=5(r), these expressions are reduced to the single-folded (SF) form. In the frame 
of the LDA the SF procedure leads to the nuclear matter approximation (NM) ||. 

In this work we have evaluated Eqs. (1) in momentum space Jl|,|^Q] by means of the 
detailed procedure developed by Khoa et al. @, and the following approximations. The 
effective ^^-interaction has been obtained by using the isoscalar and isovector components 
of the direct Vq , and exchange parts Vq Q x , of the M3Y interaction based on the results 
of the (yf-matrix calculations using either Reid or Paris iViV-potential The energy- and 
density-dependent M3Y effective iVTV-interaction has been taken in the form 

C(oif(P, E, r) = F(p) g(E) (r). (2) 

Three forms of the effective iViV-interaction density dependence, namely DDM3Y1 and 
BDM3Y1 M 



n _( C[l + aexp(-pp)], for DDM3Y1, 
*W-\C{l-ap), for BDM3Y1, 



(3a) 
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and the third one HLM3Y [10| with different density dependence for the isoscalar and 
isovector components 



F(p) 



fo(p) = s [l + a? (p/p ) 1/3 + a°(p/p ) 2/3 + ^(p/po) 1 ] 
f T {p) = s T [l + al(p/p ) 1/3 ] 



(3b) 



where involved in the present work. The energy- dependent factor is taken to be a linear 
function of the incident energy per nucleon E of the form g(E)=l-jE. The coefficient 7 has 
the values 0.002 for the M3Y-Reid interaction, and 0.003 for the M3Y-Paris interaction. 

The frozen- density approximation [|],[8],|§ has been adopted for the overlap density which 
enters the explicit form of the density-dependent factor F(p), being taken as the sum of 
the densities of the two colliding nuclei at the midpoint of the intranucleonic separation. A 
widely used approximation for the calculation of the knock-on exchange term of the folded 



potential is that proposed by Campi and Bouyssy [16], which preserves the first term of the 
expansion given by Negele-Vautherin [17] for the realistic density-matrix expression 



p(R,R + s) 



p(R+^)ji 



k av (R 



(4) 



where ji(x) = 3 (sin x — x cos x)/x 3 , and k av defines the average relative momentum [[16 

5 / I 



3p(r) V T(r) " 4 V2p{r] 



(5) 



The last quantity is a function of the density distribution p(r), and the approximation used 
for the kinetic-energy density r(r) for each participant in the interaction. 

The Fermi- type nucleon density distribution corresponding to the parameters po, R1/2, 
and a given by either the Skyrme-Hartree-Fock systematics |2(J or Negele parametrization 
pl| are compared in Figs. 1(a) and (b). The Fermi distribution for the target nucleus 
density, obtained by using the Negele parametrization |H] is shown in Fig. 1(c). 

Different forms of the kinetic-energy density r(r) have been adopted for the light pro- 
jectiles and for the target nuclei. Thus, the modified Thomas-Fermi (MTF) approximation 
|18| has been considered for light projectiles 



r(p) = ap(r)°i + (3 



I(VP)P 
P 



(6) 



The values of the parameters a and (3 depend on the number of the constituent nucleons of 
the projectiles. Following the analysis of Krivine and Treiner f!8 |, we used for the nucleon 
as a projectile a = and (3 = 1/4. For heavier nuclei and also within the MTF frame we 
have a = 3(37r 2 ) 2 / 3 /5. Consequently, the following expression have been obtained for k av (r) 
in the case of light projectiles (here nucleons) 



MTF, 



5|Vp(r)| 2 5V 2 p(r) 



12p 2 (r) 12p(r) 



(7a) 



while for heavier nuclei we have 
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■„MTF i 



k F (rf + 



5|Vp(r)| 2 _ 5V 2 p(r) 
12p 2 (r) 12p(r) 



(7b) 



where fc_F(r) = [37r 2 p(r)/2] 1//3 is the local Fermi momentum. 

On the other hand, we have also used the extended Thomas- Fermi (ETF) approximation 



22 



T{p) = ap{r)i +/5KM£ + 7 v 2 p, 

P 



where the widely-used values of the parameters a, (3, and 7 [22 
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5 V " ' 1 r 36' 
lead to the following expression for the average relative momentum 

btf,^_\l ^2 , 5| Vp(r) | 2 , 5V 2 p(r) 



(8) 



(9) 



k 



k F (rY + 



108p 2 (r) 36p(r) 



(10) 



The MTF kinetic-energy densities for neutrons and protons are shown in Figs. 1(d) and 
(e), respectively. The kinetic-energy densities derived in the ETF and MTF approximations 
are in close agreement for 93 Nb, as seen in Fig. 1(f). These results prove the consistency of 
the parametrizations involved. The related k av (r) curves are shown in Figs. 1(g), (h) and (i) 
for nucleons and 93 Nb, respectively. The increased nucleon fc a „-values at smaller radii follow 
the behaviour of the V 2 p(r) for a Fermi distribution, whose importance in the nucleon case 
results from Eq. fl7a|). 

Furthermore, the local incident momenta ki(r) calculated for 20 MeV neutrons and 
protons on 93 Nb by using the phenomenological OMP of Lagrange-Lejeune [19| as well 



as the actual DF potential are also shown in Figs. 1(g) and (h). The minor differences in 
both shape and magnitude support the use of the Negele-Vautherin approximation fl7j as 
Eq. dj) and the MTF approximation for the calculation of k av {r). The k^ (r) curve for 
the target nucleus 93 Nb decreases at larger r, however slower than k F {r). 

The real OMPs for incident neutrons and protons on 93 Nb, calculated by means of (i) the 
single-folding procedure within the LDA, (ii) the complete SF procedure, and (iii) the full 



DF, are compared with both the phenomenological OMP ]T9| and the JLM parametrization 
|| in Fig. 2. The DDM3Y1-Paris effective iViV-interaction, the Fermi density distribution 
p0| and the MTF approximation for k av (r) have been used in this respect. Moreover, the 
corresponding volume integrals Jy and rms radii (r 2 ) 1//2 are given in Table 1. The real 
part of the JLM potential has been chosen for the present analysis since the parameters of 
the DDM3Y1 interaction [|8]|| was obtained by fitting of the JLM results for the nucleon 
scattering 0|||§. On the other hand, the well-known OMP parameter set of Lagrange- 
Lejeune |0| has also been computed for both kinds of projectiles on 93 Nb, close to the JLM 
analysis. 

First, Fig. 2 shows that more diffuse potentials with larger (r 2 ) 1 / 2 values are obtained 
for neutrons as well as for protons when the folding procedure increases in accuracy going 
from NM to DF cases. Actually, the comparison in Figs. 2 (a),(b) and (d),(e) of the real 
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OMP obtained by using the NM and SF methods evidences the effects of the target nuclear 
density and the corresponding average relative momentum. Next, the comparison with Figs. 
2(c) and (f) shows even stronger effect of the nucleon density distribution and the average 
relative momentum. On the other hand, the different choices of the target k av (MTF or 
ETF) have practically no effect on the single-folded real OMP, the corresponding curves in 
Figs. 2(b) and (e) being overlapped. 

Second, the comparison between the real parts of the OMPs shown in Fig. 2 to the 
phenomenological values |0| and the JLM parametrization [|| is completed by considering 
the energy dependence of the corresponding volume integrals. From Fig. 3 it is seen that the 
phenomenological and both calculated Jy values from the SF and the DF approach agree 
rather well. In view of the increasing agreement for the depths (Fig. 2) and the (r 2 ) 1 / 2 
values (Table I) in going from the SF to DF procedures these comparisons can be considered 
an important test of the present calculations. A comment concerning the proton case is 
necessary. There, the comparison to the JLM parametrization is less significant, because 
the latter accounts only for the isoscalar component of the full real OMP. 

Third, we have considered in the present work three different density-dependent versions 
DDM3Y1, BDM3Y1 and HLM3Y of the M3Y-Paris and M3Y-Reid iViV-interactions. The 
double-folded real OMPs corresponding to the former interaction for DDM3Y1 and BDM3Y1 
are shown in Fig. 4 for the neutron as well as the proton cases. It results that the two density 
dependences lead to almost identical folded real OMPs provided that a Fermi-distribution 
of the nucleon density is used. Also, the calculated real part of the OMP using HLM3Y 
interaction (JTUJ] lead to comparable results with those obtained by DDM3Y1 and BDM3Y1. 
The radial dependence of the real OMP corresponding to the Paris-HLM3Y effective NN- 
interaction, however evidences enhanced values in the surface region as can be seen from the 
rms radii values in Table 1. On the other hand, an obvious difference exists between the 
DF real potentials corresponding to the simple M3Y-interaction and the density-dependent 
versions, respectively. 

The density dependence of the effective iViV-interaction improves the agreement of the 
folded potentials with the phenomenological ones by reducing the DF potential values at 
small values of r where the highest overlap of the interacting projectile and target-nucleus 
densities takes place. However, the simple M3Y folded potentials are less diffuse in the 
nuclear surface region, e.g. they have smaller (r 2 ) 1 / 2 values than the real OMPs obtained 
by means of the energy- dependent iViV-interactions (see Table 1). 

The importance of the realistic effective ^^-interactions is pointed out by showing also 
in Fig. 4 the DF real potentials calculated by means of the usual 1 fm range Yukawa 
interaction with the V strength values obtained by Watanabe et al. for these systems p| . 
i.e., V =40 MeV for n+ 93 Nb and 46.8 MeV for p+ 93 Nb. On the other hand, there are 
also shown the results obtained by using the equivalent 1 fm range Yukawa interaction 
with Vq strengths corresponding to the volume integral values provided by the Paris-M3Y 
interaction. Thus, the improvement in comparison with the 1 fm range Yukawa interaction 
obtained just using the simpler M3Y interaction is obvious. Rather similar effect has the use 
of the equivalent Vq values, corresponding to the Paris-M3Y volume integrals. Obviously, 
the supplementary density dependence increases the agreement with the phenomenological 
optical potentials [see Figs. 2(c) (f)] considerably. 

Finally, one comment may concern the effect of the ^^-interaction type, i.e. of the Paris 
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and Reid potentials. As seen from Table 1, the Jy and (r 2 ) 1 / 2 values for the corresponding 
neutron and proton folded potentials are close to each other, independent of using the NM, 
SF, and DF. Thus one may conclude that both NN interactions lead to reliable results in 
folding calculations. 

In conclusion, this analysis has shown that the double folding procedure involving the 
nucleon-density Fermi distribution and average relative momentum leads to real OMPs at 
lower energies in good agreement with the phenomenological OMPs as well as with the 
JLM parametrization for both neutron- and proton-nucleus scattering. The analysis of 
the average effective iViV-interaction strength by using the actual folded real OMP and 
following [0] is in progress. Since this quantity is still considered as the only free parameter 
of the quantum-statistical studies of the multistep reactions, the corresponding search for 
the realistic effective iViV-interaction at low energies is thus of further interest. 

The comparison of the full folding calculations to the results obtained with the strongly 
simplified single Yukawa description is very instructive showing clearly the restrictions of a 
"simple" parameterization. Apparently, the depth of the potential is overestimated because 
of the lack of density dependence, as one should expect. For inelastic transitions at a 
fixed incident energy one might compensate that shortcoming by simply rescaling the Vo 
parameter such that the "local" interaction strength in the reaction region is accounted 
for in the average. However, the penetration and therefore the interaction region changes 
with incident energy and Q-value. Very likely, the observed energy dependence of V is 
actually produced to a large extent by the underlying density dependence which was taken 
into account explicitly in this work. 

A clearer account of that effect could be obtained from elastic and inelastic scattering 
reaction calculations (OMP and DWBA) at different incident energies where one tries to 
reproduce the cross sections derived with density dependent interactions by OMPs and form 
factors from a single Yukawa interaction. 
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Table 1. Volume integrals (in MeV-fm 3 ) and radii (r 2 ) 1 / 2 (in fm) for the real part 
of folded OMP calculated by means of the nuclear matter (NM), single-folding (SF) and 
double-folding (DF) procedures, phenomenological OMP [pl| and JLM parametrization ||. 
Interaction NM SF DF Phenomenological 

Jv {r 2 ) 1 ' 2 Jv (r 2 ) 1/2 Jv (r 2 ) 1/2 Jv (r 2 ) 1 / 2 



neutrons 



JLM 404.5 4.652 



398.6 4.930 



Paris-M3Y 383.0 4.399 402.0 4.692 398.7 5.057 

Paris-HLM3Y 448.3 4.812 470.3 5.072 422.4 5.270 

Paris-HLM3Y* 452.9 4.857 

Paris-BDM3Y1 385.7 4.527 404.9 4.808 384.9 5.168 

Paris-DDM3Y1 395.7 4.567 415.3 4.846 392.2 5.194 

Reid-DDM3Y1 360.2 4.568 383.7 4.854 367.0 5.206 



protons 

Paris-M3Y 429.3 4.413 452.7 4.692 451.4 5.075 

Paris-HLM3Y 501.1 4.793 528.2 5.046 483.0 5.268 
Paris-HLM3Y* 507.7 4.839 

Paris-BDM3Y1 433.1 4.539 456.9 4.808 438.0 5.183 

Paris-DDM3Y1 444.2 4.579 468.6 4.846 446.3 5.209 

Reid -DDM3Y1 422.0 4.584 453.0 4.865 436.4 5.234 



470.1 4.930 



*The Hartree-Fock nuclear density for 93 Nb has been used in this case [[2(J 
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FIGURE CAPTIONS 



FIG. 1. Nuclear density distributions, kinetic-energy densities r(r) (MTF, solid curves) and 
(ETF, dotted curve), the local incident momenta ki corresponding to the 20 MeV 
incident energy, and MTF average relative momenta k av (r), for (a,d,g) neutron, (b,e,h) 
proton, and (c,f,i) the target nucleus 93 Nb, respectively. The Fermi-distribution as 
given by Lenske |2(J (solid curves) or Negele parametrization (dotted curves) 
are used for nucleon density distribution, and only the latter for the target nucleus. 
The ki values correspond to the phenomenological OMP of Lagrange-Lejeune (LL) 
[[0| (dashed-dotted curves) and to the DF potential (dashed curves), while the MTF 
average relative momentum (solid curve), the ETF k av (solid curve), and the local 
Fermi momentm (dotted curves) are shown for the taget nucleus. 

FIG. 2. The radial dependence of the real OMP calculated with the Paris-DDM3Y1 effective 
iVW-interaction in the MTF approximation, by using the nuclear matter approxima- 
tion (NM), single- folding (SF) and double-folding (DF) procedures (solid curves), com- 
pared with the phenomenological OMP |1J (dashed curves) and JLM parametrization 
[||] (dotted curves), for (a,b,c) neutrons and (d,e,f) protons of 20 MeV incident energy 
on the target nucleus 93 Nb. 

FIG. 3. The incident-energy dependence of the volume integrals of real OMP calculated with 
the Paris-DDM3Y1 effective iViV-interaction in the MTF approximation, by using 
the nuclear matter approximation (NM), single-folding (SF) and double- folding (DF) 
procedures (solid curves), as well as with the phenomenological OMP |19| (dashed 
curves) and JLM parametrization || (dotted curves), for (a,b,c) neutrons and (d,e,f) 
protons incident on the target nucleus 93 Nb. 

FIG. 4. The same as Fig. 2, for (a) neutrons and (b) protons, and corresponding also to the 
effective NN interactions M3Y (dotted curves), BDM3Y1 (dashed curves), 1 fm range 
Yukawa with the Vo values either obtained through FKK analysis [^3|] (dotted-dashed 
curves) or equivalent (double dotted-dashed curves) to the M3Y interaction (see text). 
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